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Abstract 

In this paper a detailed simulation of irradiated pixel sensors was used to investigate the effects of radiation damage 
on the position determination and optimize the hit reconstruction algorithms. The simulation implements a model 
of radiation damage by including two defect levels with opposite charge states and trapping of charge carriers. 
The simulation shows that a position resolution below 15 /J,m along the CMS r — <j> plane can be achieved after an 
irradiation fluence of 5.9 x 10 14 n eq /cm 2 . In addition, we show that systematic errors in the position determination 
can be largely reduced by applying rj corrections. 



1. Introduction 

The CMS experiment, currently under construc- 
tion at the Large Hadron Collider (LHC) will in- 
clude a silicon pixel detector [1] to allow tracking 
in the region closest to the interaction point. The 
detector will be a key component for reconstruct- 
ing interaction vertices and heavy quark decays in 
a particulary harsh environment, characterized by 
a high track multiplicity and heavy irradiation. At 
the full LHC luminosity the innermost layer, with 
a radius of 4.3 cm, will be exposed to a particle 
fluence of 3 x 10 14 n eq /cm 2 /yr. 

In order to evaluate the effects of irradiation and 
optimize the algorithms for the position determi- 
nation a detailed simulation of the pixel sensors 
was implemented. In [2] we have proved that it is 
possible to adequately describe the charge collec- 



tion characteristics of a heavily irradiated silicon 
detector in terms of a tuned double junction model 
which produces a doubly peaked electric field pro- 
file across the sensor. The modeling is supported 
by the evidence of doubly peaked electric fields 
obtained from beam test measurements and pre- 
sented in [3] . The dependence of the modeled trap 
concentrations upon fluence was presented in [4] 
and the temperature dependence of the model was 
discussed in [5]. In this paper the simulation was 
used to study the position determination in irradi- 
ated pixel sensors. 

This paper is organized as follows: the sensor 
simulation is described in Section 2, in Section 3 the 
determination of the hit position in irradiated pixel 
sensors is discussed. The results are presented in 
Section 4 and the conclusions are given is Section 5. 
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2. Sensor simulation 

The results presented in this paper rely upon a 
detailed sensor simulation that includes the mod- 
cling of irradiation effects in silicon. The simula- 
tion, pixelav [2,6,7], incorporates the following 
elements: an accurate model of charge deposition 
by primary hadronic tracks (in particular to model 
delta rays) ; a realistic 3-D intra-pixel electric field 
map; an established model of charge drift physics 
including mobilities, Hall Effect, and 3-D diffusion; 
a simulation of charge trapping and the signal in- 
duced from trapped charge; and a simulation of 
electronic noise, response, and threshold effects. 
The intra-pixel electric field map was generated 
using TCAD 9.0 [8] to simultaneously solve Pois- 
son's Equation, the carrier continuity equations, 
and various charge transport models. 

The simulated devices correspond to the base- 
line sensor design for the CMS barrel pixel detec- 
tor. The sensors are "n-in-n" devices, designed to 
collect charge from n+ structures implanted into 
n- bulk silicon. The simulated samples were 22x32 
arrays of 100x150 fim 2 pixels. The substrate was 
285 /im thick, n-doped silicon. The donor concen- 
tration was set to 1.2 x 10 12 cm~ 3 corresponding 
to a depletion voltage of about 75 V for an unir- 
radiated device. The 4 T magnetic field was set as 
in the CMS configuration and the sensor temper- 
ature to -10° C. The simulation did not include 
the "punch-through" structure on the n + implants 
which is used to provide a high resistance connec- 
tion to ground and to provide the possibility of on- 
wafer IV measurements. 

The effect of irradiation was implemented in the 
TCAD simulation by including two defect levels in 
the forbidden silicon bandgap with opposite charge 
states and trapping of charge carriers. The activa- 
tion energies of the donor and acceptor traps were 
set to (E v + 0.48) eV and (E c - 0.525) eV, re- 
spectively, where Ey and Ec are the valence and 
conduction band energy level, respectively. The 
trap densities and the capture cross sections for 
electrons and holes were obtained by fitting the 
model to beam test data as described in [2,4]. The 
simulated irradiation fluences were <& = 2 x 10 14 
n cq /cm 2 and $ = 5.9 x 10 14 n oq /cm 2 and the re- 



verse bias was set to 200 V and 300 V, respectively. 

The electric field profile as function of the sen- 
sor depth is shown in Fig. 1 (a) . The field has max- 
ima at the detector implants and a minimum near 
the midplane which shifts towards the p + implant 
at lower fluences. The dependence of the space 
charge density upon the detector depth is shown 
in Fig. 1(b). Before irradiation the sensor is char- 
acterized by a constant and positive space charge 
density across the sensor bulk. At a fluence $ = 
2 x 10 14 n eq /cm 2 the device shows a negative space 
charge of about —6 x 10 13 cm~ 3 for most of its 
thickness, a compensated region corresponding to 
the electric field minimum and a positive space 
charge density close to the backplane. The space 
charge density and electric field near the p + im- 
plant increase with the irradiation fluence. 
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Fig. 1. The z-component of the simulated electric field (a) 
and space charge density (b), resulting from the two-trap 
model, arc shown as a function of z for unirradiated and 
irradiated devices. 



3. Hit reconstruction 

The spatial resolution of the pixel sensors is 
mainly determined by the readout pitch and charge 
sharing between neighboring cells. Pixels have a 
weak capacitive coupling and charge sharing is 
mainly due to diffusion and drift of charge carriers 
under the combined effect of the magnetic and 
electric fields. After irradiation, free carriers trap- 
ping produces an inhomogeneous charge collection 
across the bulk and charge sharing between neigh- 
boring pixels becomes strongly non linear on the 
impact position. In addition, the beneficial effect 
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of the Lorentz deflection is reduced when a higher 
bias voltage is applied to provide a sufficient drift 
field. In what follows we discuss measurements of 
the sensor spatial resolution along the r — <fi di- 
rection, where the charge drift is affected by the 
Lorentz deflection. 

To reconstruct hits pixels with charge above 
2000 electrons were selected and clusters were 
formed by adjacent pixels above threshold. Both 
side and corner adjacent pixels were included in 
the cluster. In addition, clusters adjacent to the 
matrix border were excluded. Figure 2 shows the 
definition of the track impact angle a with respect 
to the sensor plane along the x direction. The 
track is orthogonal to the sensor plane along the 
axis orthogonal to x. The magnetic field produces 
a Lorentz shift L = Ttan(Oi) towards the right 
direction. Thus, the total charge width is given by 



W =i-Ttan(a). 



(1) 



The cluster is projected along the x axis by sum- 
ming the charge collected in the pixels with the 
same x coordinate. If the cluster is given by a single 
pixel its position is given by the pixel center. For 
larger clusters the hit position is calculated with a 
center-of-gravity algorithm. 
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Fig. 2. Determination of the impact position in the trans- 
verse plane. 

To further improve the spatial resolution for 
events in which charge is shared among several 
pixels the so-called 77-correction is introduced [9]. 
As we will discuss, the correction is particularly 
effective on irradiated devices, where the effects 
of inhomogeneous charge collection are larger. As- 
suming that the number of particles crossing the 
sensor is uniformely distributed along x we should 
expect the reconstructed position within a pixel 
to be uniformely distributed. We define 77 as the 
non-integer part of the reconstructed position. 



Figure 3 shows the distribution of 77 for all events, 
where 77 = corresponds to the center of the pixel 
cell and 77 = ±0.5 to the borders. The measured 
distribution is almost flat in the pixel regions 
closer to the pixel borders and it shows a dip at 
the center. The peak at r\ = is due to single pixel 
clusters. For each 77 we associate a corrected value 
given by the function 



J\ 5 dN/dfjdfj l 
f° 5 5 dN/dfj dfj 2 



(2) 



where 77 is in pixel units. The ^(77) function is 
shown in Fig. 3. The corrected position is calcu- 
lated by adding F{q) to the integer part of the re- 
constructed position. 



4. Cluster size and position resolution 

Figure 4 shows the fraction of events for different 
cluster sizes as function of the impact angle a and 
for different irradiation fluences. At negative angles 
the Lorentz shift and the geometrical sharing term 
sum up giving large W values. For perpendicular 
tracks (a = 0°) the total charge width is given by 
the Lorentz shift and then decreases for a > 0°. 
Clusters are larger than one pixel in the majority 
of cases, however at high irradiation fluences the 
higher bias voltage produces a narrower Lorentz 
deflection and, consequently, a smaller cluster size. 

The combination of carrier trapping and higher 
bias voltage produces smaller cluster sizes after 
heavy irradiation. The average cluster size for an 
unirradiated detector with Vbi as =200 V and per- 
pendicular tracks is 2.0. After a fluence of $ = 
5.9 x 10 14 n cq /cm 2 and Vf,i as =300 V the value is 
1.8. 

The position resolution was estimated by com- 
paring the hit position, x rec , reconstructed as de- 
scribed in Section 3 with the true impact position, 
Xtruei generated by the simulation. The residuals 
were defined as x res = x rec — Xtrue and the position 
resolution was given by the width of a Gaussian fit 
of the residual distribution. The position resolution 
for tracks perpendicular to the sensor is summa- 
rized in Table 1 for clusters of two pixels. The posi- 
tion resolution for all events is shown in the last col- 
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(a) (b) 
Fig. 3. (a) Distribution of the reconstructed impact position within a single pixel for perpendicular tracks, (b) Correction 
function F(r/) (see Eq. 2). 




Fig. 4. Event fractions for different cluster sizes as function of the impact angle a for an unirradiated sensor (a) and for 
sensors irradiated to $ = 2 X 10 14 n cq /cm 2 (b) and <3? = 5.9 X 10 14 n eq /cm 2 (c). 
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Table 1 

Position resolution for a = 0° at different irradiation fluences and bias voltages. The third and fourth column show the 
position resolution of the two-pixel clusters without and with r) correction, respectively. The last column shows the resolution 
for all events, where the position of the two-pixel clusters was 77 corrected. 



umn, where the position of the two-pixels clusters 
was rj corrected. The simulation shows that posi- 
tion resolution below 15 can be achieved even 
after heavy irradiation. In addition, the precision 
can be improved by correcting the reconstructed 
position as described in Section 3. 

The size of the rj corrections become particu- 
larly significant after irradiation and for narrow 



charge widths. Figure 5 shows the residuals distri- 
bution for clusters of two pixels, simulated for a 
sensor irradiated to $ = 5.9 x 10 14 n oq /cm 2 and 
for tracks with a = 20°. The distribution before 
correction (Fig. 5(a)) is not described by a single 
Gaussian and is affected by large systematic errors 
which depend on the interpixel hit position. The 
distribution in Fig. 5(b) shows that the systematic 
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Fig. 5. Residuals distribution for tracks with a = 20° and a sensor irradiated to $ = 5.9 X 10 14 n oq /cm 2 . The distributions 
are calculated without (a) and with ^-corrections (b) and are not corrected for the Lorentz shift due to the magnetic field. 
The simulated data points are represented by the markers and the continuous line is a double-Gaussian fit to the distribution. 



errors can be largely reduced by applying the rj- 
correction. 



5. Conclusions 

In this paper a detailed simulation of the sili- 
con pixel sensors for the CMS tracker was used 
to estimate the effects of radiation damage on the 
position resolution. The simulation, incorporating 
a double junction model of radiation damage and 
trapping of charge carriers, provides a good de- 
scription of the charge collection measurements 
in the fluence range from 0.5 x 10 14 n eq /cm 2 to 
5.9 x 10 14 n cq /cm 2 . 

The simulation shows that a position resolution 
below 15 fim can be achieved for perpendicular 
tracks in the CMS trasverse plane even after heavy 
irradiation. In addition, we show that the position 
resolution can be improved by applying r\ correc- 
tions. 
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